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ABSTRACT 

We measure the far-infrared emission of the general quasar (QSO) population using Planck ohservations of the Baryon Oscillation 
Spectroscopic Survey QSO sample. By applying multi-component matched multi-filters to the seven highest Planck frequencies, we 
extract the amplitudes of dust, synchrotron, and thermal Sunyaev-Zeldovich (SZ) signals for nearly 300 000 QSOs over the redshift 
range 0.1 < z < 5. We hin these individual low signal-to-noise measurements to obtain the mean emission properties of the QSO 
population as a function of redshift. The emission is dominated by dust at all redshifts, with a peak at z ~ 2, the same location 
as the peak in the general cosmic star formation rate. Restricting analysis to radio-loud QSOs, we find synchrotron emission with 
a monochromatic luminosity at 100 GHz (rest-frame) rising from Lsynch = 0 to 0.2LoHz“' between z = 0 and 3. The radio-quiet 
subsample does not show any synchrotron emission, but we detect thermal SZ between z = 2.5 and 4; no significant SZ emission is 
seen at lower redshifts. Depending on the supposed mass for the halos hosting the QSOs, this may or may not leave room for heating 
of the halo gas by feedback from the QSO. 

Key words, cosmology: observations - large-scale structure of Universe - quasars: general - galaxies: clusters: general - methods: 
data analysis - methods: statistical 


1. Introduction 

Quasars (or quasi stellar objects, QSOs) occupy a special place 
in large-scale structure and galaxy evolution (Kormendy & Rich- 
stone 1995). They are among the most luminous extragalactic 
sources, and, as such, have become the focus of many cosmo¬ 
logical surveys such as the 2dF Quasar Redshift Survey (2QZ; 
Groom et al. 2001) or the successive iterations of the Sloan Dig¬ 
ital Sky Survey (SDSS TIV; York et al. 2000; Eisenstein et al. 
2011), where they provide a unique insight into the formation of 
structure on large scales and at high redshift (White et al. 2012). 
Quasars are likely powered (Salpeter 1964; Lynden-Bell 1969) 
by accretion of nearby matter onto supermassive black holes 
(SMBH, Rees 1984), whose evolution appears closely linked 
to the general cosmic star formation rate (Madau & Dickinson 
2014), in particular in galaxies that contain a massive bulge (and 
therefore a massive central black hole) and a gas reservoir (Nan- 
dra et al. 2007; Silverman et al. 2008). The link is a clue to 
galaxy formation, its significance emphasized by the observed 
relation between the SMBH mass and galaxy properties (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000). In addition, galaxy 
formation models must evoke strong feedback from AGN to ex¬ 
plain the observed properties of massive galaxies and to avoid 
the overcooling catastrophe (Croton et al. 2006; Bower et al. 
2006; McNamara & Nulsen 2007; Somerville et al. 2008; Fabian 
2012; Blanchard et al. 1992) 

Quasar environments give us a look at these powerful en¬ 
gines, and millimeter/submillimeter observations can offer a par¬ 
ticularly revealing view of their effects: at these frequencies it is 


possible to study the cooler dust emission associated with star 
formation in the host, look for synchrotron emission from ener¬ 
getic particles, and, perhaps most pertinently, to measure energy 
feedback through observation of the thermal Sunyaev-Zeldovich 
(tSZ, Sunyaev & Zeldovich 1970, 1972) effect, a direct probe of 
the thermal energy contained in surrounding gas. Samples ob¬ 
served in this waveband with large ground-based facilities typ¬ 
ically consist of, at most, several tens of objects (Omont et al. 
2001; Isaak et al. 2002; Omont et al. 2003; Priddey et al. 2003; 
Beelen et al. 2006; Wang et al. 2007; Omont et al. 2013). Oper¬ 
ating at these same frequencies, cosmic microwave background 
(CMB) surveys cover large sky areas and present the opportunity 
of studying much larger samples, albeit with much less detail 
of individual objects. We can, however, determine the average 
properties of large representative populations, a valuable com¬ 
pliment to the more detailed examinations. 

The Wilkinson Microwave Anisotropy Probe (WMAP, Ben¬ 
nett et al. 2003) and Planck missions (Tauber et al. 2010; Planck 
Collaboration et al. 2011a) are ideal for this purpose because 
of their all-sky coverage, including the entire Sloan Digital Sky 
Survey (SDSS, York et al. 2000) area. Cross-correlating WMAP 
data with photometric QSOs from SDSS Data Release 3 (DR3), 
Chatteijee et al. (2010) finds evidence at 2.5cr for the tSZ effect 
from QSO environments. Ruan et al. (2015) improved the signif¬ 
icance using a publicly available tSZ map (Hill & Spergel 2014) 
constructed from the Planck mission dataset and 26 686 spec¬ 
troscopic QSOs from SDSS DR7. They concluded that the total 
thermal energy feedback into surrounding gas was significantly 
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larger than expected according to galaxy formation models, 
which is consistent with the findings of Chatterjee et al. (2010). 
Planck’s wide frequency coverage, from 30 GHz to 857 GHz, 
is a distinct advantage for disentangling the different sources 
of emission from the QSO environment, and was implicitly ex¬ 
ploited by Ruan et al. (2015) when using the tSZ map published 
by Hill & Spergel (2014). 

Although not all-sky, the Atacama Cosmology Telescope 
(ACT) surveyed several hundred square degrees near the celes¬ 
tial equator, overlapping AGN and QSO samples in the North. 
Gralla et al. (2014) stacked millimeter maps from ACT over this 
area of radio-loud AGN seen in the FIRST and NVSS radio sur¬ 
veys to find a tSZ signal at 5cr significance. More recently, and 
in parallel to our present study, Crichton et al. (2015) combined 
ACT and Herschel data of radio-quiet quasars in the SDSS DR7 
and DR 10 samples, detecting the tSZ signal at 3 - 4cr signih- 
cance. 

In this paper, we use the full-mission Planck frequency maps 
to examine the emission properties of the SDSS-III Baryon 
Acoustic Oscillation (BOSS, Dawson et al. 2013) DR12 QSO 
sample. We extract not only the tSZ signal from the population, 
but also dust and synchrotron emission by stacking measure¬ 
ments made with a set of matched biters directly applied to the 
seven highest frequency Planck channels. The joint extraction 
enables us to study not only the gas thermal environment, but 
star formation conditions and the production of energetic parti¬ 
cles. Moreover, it improves analysis robustness by giving infor¬ 
mation on correlations between the observed signals, compared 
to the use of a tSZ map. Ruan et al. (2015) were in fact faced 
with the difbculty of correcting the SZ map for contamination 
from dust emission. 

We extend our multi-matched biter (MMF, Melin et al. 2006, 
see also Herranz et al. (2002)) formalism to simultaneously ex¬ 
tract two or three emission components, following development 
started in Planck Collaboration et al. (2013). Stacking the vari¬ 
ous signals to study their evolution with redshift, we reconstruct 
for the brst time a picture of the evolution of the dust signal, 
the synchrotron emission, and the tSZ effect between z ~ 0 and 
z ~ 5 for the general QSO population. 

In Sect. 2, we describe the BOSS and Planck datasets. We 
detail the analysis methods and tools in Sect. 3. In Sect. 4, we 
apply them on simulations. Results on Planck data are given 
in Sect. 5. We discuss our results in Sect. 6 and conclude in 
Sect. 7. Throughout, we use the spatially bat base ACDM cos¬ 
mology from Planck Collaboration et al. (2015d): Hq - 67.27 = 
blOOkm/s/Mpc, = 0.3156, Q.bh^ - 0.02225 and erg — 0.831. 

2. Data 

2.1. BOSS quasars 

One of the major goal of the Sloan Digital Sky Survey-III 
(SDSS-III) Baryon Oscillation Spectroscopic Survey (BOSS) is 
the production of a QSO catalogue to detect the BAO scale in the 
Lyman-a forests at redshift z~ 2.5. A brst detection was made 
on the DR9 QSO catalogue (Pms et al. 2012) by Busca et al. 
(2013), conbrmed latter by Delubac et al. (2015). In this study, 
we use the recently published DR 12 QSO catalogue' (Alam et al. 
2015). 

The sources are detected and selected with the CCD imaging 
Camera installed in the Sloan Foundation 2.5 m Telescope 

' http: //WWW. sdss.org/drl2/algorithins/ 
boss-dr12-quasar-catalog/ 



Fig. 1. Redshift distribution of the BOSS DR12 QSO sample. The dis¬ 
tribution presents two peaks, at z ~ 0.8 and z ~ 2.3, due to the QSO 
target selection process (see text). The regular binning used throughout 
this paper (Az = 0.5) is shown as the vertical dotted lines. 

at Apache Point observatory. New Mexico. The spectra of 
each source is measured with the BOSS spectrograph which 
covers wavelength between 3600 A and 10 000 A. Spectra 
measurement leads to the computation of the spectroscopic 
redshift and the checking of the nature of the source. 

The DR12 QSO catalogue contains 297 301 objects. We re¬ 
move 5880 QSO falling outside the Planck 65% mask, in a re¬ 
gion strongly contaminated by the Milky Way dust. Finally, we 
also reject 256 QSO because they are at low redshift z < 0.1 
so may be partially resolved by Planck , because they belong to 
the high redshift (z>5) population or because they have a bad 
estimate of the magnitude in g band (g < 0). We thus use a cat¬ 
alogue of 291 165 QSO in the redshift range 0.1 < z < 5. The 
distribution of the QSO in term of redshift is not bat and reach 
two notable maxima at z=0.8 and z=2.3 (Fig. 1). 

The target selection was developed to select quasars with an 
observable Ly-a forest (i.e quasars with z>2.2). However, de¬ 
generacies in the color-redshift relation of quasars led to the 
selection of low-z quasars in BOSS (Fig. 1). The quasars at 
z ~ 0.8 have Mgll 42800 line at the same wavelength as Ly-a 
at redshift z ~ 3.1, giving these objects similar broadband col¬ 
ors, while the large number of objects at z ~ 1.6 is due to the 
confusion between /11549 C-IV line and Ly-a at z ~ 2.3 (Ross 
et al. 2012). In conbast, the selection based on the intrinsic vari¬ 
ability of quasars gives a more uniform distribution in redshift 
(Palanque-Delabrouille et al. 201 1). 

2.2. Planck maps 

The Planck satellite was launched in May 2009 (Tauber et al. 
2010). After traveling to the Earth-Sun Lagrange point L2, it 
scanned the entire sky continuously from August 2009 to Octo¬ 
ber 2013 in nine frequency bands ranging from 30 to 857 GHz. 
The primary goal of the mission was to study the primary 
CMB anisotropies, but its large frequency coverage also en¬ 
ables unique Galactic and extra-galactic astrophysical studies. In 
particular, the High Frequency Instrument (HFI, Lamarre et al. 
2010; Planck HFI Core Team et al. 201 1), covering bands from 
100 to 857 GHz, is ideally suited for Sunyaev-Zeldovich science. 
Moreover, the three highest frequencies (353 to 857 GHz) of the 
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instrument pick up Galactic and extra-galactic dust emission. 
The Low Frequency Instrument (LFI, Bersanelli et al. 2010; 
Mennella et al. 2011), with channels at 30, 44, 70 GHz, is sensi¬ 
tive to synchrotron, free-free, and spinning dust emission. 

We use the seven highest frequency (70, 100, 143, 217, 
353, 545, and 857 GHz) full-sky maps from the 2015 data re¬ 
lease (Planck Collaboration et al. 2015a,b,f). These bands are 
the best suited to study the dust, gas, and synchrotron emission 
from high redshift QSO environments. The substantially larger 
beams of the remaining 30 and 40 GHz channels make them less 
well adapted for studying what are essentially point sources. We 
divide each Planck HEALPix^ all-sky map into 504 overlapping 
flat tangential patches of 10 x 10 degrees with a 1.72 arcmin 
pixel scale in order to apply our extraction algorithms described 
in Sec. 3.2. 


3. Analysis 

Determining the physical properties of the QSO environment 
from the Planck maps faces two major challenges. The first is 
the faint QSO flux in the Planck maps, far below the noise level; 
individual detection is not possible. We thus have to use a statis¬ 
tical approach similar to the one developed in Melin et al. (2011) 
and subsequently in Planck Collaboration et al. (201 ld,e, 2013). 
Details of the approach are given in Sect. 3.4.1. 

The second challenge arises from the superposition of differ¬ 
ent emission sources from the QSO in the Planck beams: emis¬ 
sion from dust located inside the host galaxy or possibly at larger 
scale in the host halo of the QSO; synchrotron emission from the 
host galaxy or from relativistic outflows of the central AGN; and 
the tSZ effect due to the hot gas surTounding the QSO and within 
the host halo. Planck does not have the spatial resolution to sep¬ 
arate these components, but its good spectral coverage enables 
us to disentangle these sources of emission. 

We separate the signals using multi-component matched 
multi-filters (MMF), an extension of the approach used in Planck 
Collaboration et al. (2013). The detailed description is given in 
Sec. 3.2.2. We will see that the QSO signal in Planck is domi¬ 
nated by dust emission, but that we also detect both synchrotron 
and tSZ signals. In practice, we proceed as follows: 

1. Assume that the QSO signal is a mixture of one, two or three 
components; 

2. Apply the adapted multi-component MMF at each QSO po¬ 
sition to obtain the amplitude of each component; 

3. Bin average (as a function of redshift or magnitude) these 
amplitudes over the QSO sample; 

4. Evaluate the ability of our model to describe the data using a 
;^^^-test (see Sect. 3.4.2). 


position X on the sky is written as 

m(x) = ^ fli G(z, x) -H n(x), ( 1 ) 

A 

where aA is the amplitude of the A component, t.i(z, x) is the as¬ 
sociated, normalized emission vector, z the redshift of the QSO, 
and n(x) is the astrophysical and instrumental noise. We center 
the QSO in the map to simplify the expressions. The compo¬ 
nent of each emission vector is the signal profile, ta (normalized 
to unity at the origin), convolved by the Planck beam, bi (nor¬ 
malized to unity at the origin), at frequency v,, and then scaled 
by the expected frequency dependance (5^), (z): 


(G); {Z, x) = (5 A)i (Q [bi * T^](x). 


( 2 ) 


We describe the dust emission frequency dependence by a 
modified blackbody, the synchrotron emission by a power law, 
and we use the non-relativistic calculation for the tSZ (e.g., 
Birkinshaw 1999; Carlstrom et al. 2002): 


(5d),(z) = 


Vi Bv,[rd/(l+z)] 


\ 857GHz / 


i 857GHz; B857GHz[7’d/(l+z)]’ 

GHzr- 

L / 






{S sz)i - By^Tcmb) 


e '- - 1 


tanh(.r:/ 2 ) 


(3) 

(4) 

(5) 


where x, = hvilkTcmh with Tcmb the temperature of the CMB 
today, By{T) represents the Planck spectrum, and /3d are the 
dust temperature and emissivity index, and as = 0.7 for the syn¬ 
chrotron index. It is the typical value for radio galaxies (Condon 
1992; Peterson 1997). We note that the first two expressions are 
unitless, so that adust and asynch carry units of brightness; on the 
other hand, Stsz carries brightness units while atsz is unitless and 
corresponds to the central Compton-y parameter. 

We adopt the same spatial profile for all signal templates, 
Ta - T(x/6f), where r is taken from Arnaud et al. (2010) based 
on the observed cluster halos gas pressure profile. We fix Of = 
0.27 arcmin, the value corresponding to a halo of mass M 500 = 
IO'^Mq at z = 2. This value is smaller than the smallest Planck 
beam (fwhm = 4.2 arcmin at 857 GHz): the bulk of the QSO 
remain unresolved by Planck. We study the sensitivity of our 
results to this fixed size in Sect. 5.5. The profiles are cut at 6 > 

50500 - 

In the following we therefore denote this universal QSO pro¬ 
file by Tqso. We are then able to define the QSO brightness, f, as 
the sum over the different components as 



(6) 


Our notation in the subsequent sections is as follows. We use 
the letter p as index to denote individual QSOs in the catalogue 
and the indices i and j to denote observation frequencies. The 
Greek indices A and p specify the nature of the emission compo¬ 
nents, e.g., A and p have as possible values dust, synch, and tSZ 
for dust, synchrotron, and tSZ signals, respectively. 


and rewrite the map at frequency i as 

OT,(x) = fi (fqso). (x) -H n,(x), (7) 

where (fq^o), (x) = [bi * Tq^oKx). 


3.1. Model of the QSO emission 

We model the QSO emission as a sum of dust, synchrotron, and 
tSZ signals. The seven-frequency column vector map, m(x), at 

^ http://healpix.sourceforge.net 


3.2. Matched filters 

3.2.1. Single frequency matched filter 

We use single frequency matched Alters to extract individual 
QSO signals {at S /N <Si 1) at each Planck frequency. The indi¬ 
vidual signals are necessary to compute thex^ statistic described 
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in Sect. 3.4.2. We also average them to obtain the mean QSO 
spectrum in Sect. 5. 

The single frequency matched filter is a linear filter de¬ 
signed to return an unbiased signal estimate with minimal er¬ 
ror (Haehnelt & Tegmark 1996; Sanz et al. 2001). No assump¬ 
tion is made concerning the frequency dependance of the signal. 
The estimated value, fi, of the true signal, /i, is 

Jd^x'P'(x)mi(x), (8) 


with different MMF to separate dust and tSZ signals. Multi- 
component filtering has also been successfully used to jointly 
extract the tSZ and kinetic SZ (kSZ) effects in Planck Collab¬ 
oration et al. (2014c). For these filters, the recovered amplitude 
of each component is unbiased with minimal variance if the as¬ 
sumption on the number and type of the components is correct. 

We consider three multi-component filtering schemes; 1) 
dust-HtSZ, 2) dust-Hsynchrotron, and 3) dust-i-tSZ-i-synchrotron 
filters. The amplitudes oa are estimated using a linear combi¬ 
nation of the partial MMFs, 


where T'i(x) is the matched filter for frequency v,. The filter T'i 
is expressed in Fourier space as 


'Pi(k) = 


(K) 


Piiik) ’ 

with the error on the estimated signal 


cr[fi] = - (Af 


f 


d^k 


|(fqso). (k)| 


- 1/2 


(9) 


( 10 ) 


where P,i(k) is the noise power spectrum of the Planck map, 
OT,(k). Because the signal-to-noise of an individual QSO signal 
is smaller than unity, we estimate P„(k) directly as the power 
spectrum of the map. 


3.2.2. Matched Multi-Filters (MMF) 


dx 




(14) 


with the partial MMF defined in Fourier space by 

cl>^(k) = p-'(k)-Vz,k), (15) 

and the matrix D computed as 


Da,ia 



d^k ty(z,k)-p-'(k)-t^(z,k) 


(16) 


The estimated amplitudes are combined together in a N- 
component vector {N being the number of signal components 
considered), a. They are correlated with covariance matrix 


= ['>■'1,^ (11) 

In summary, we use the single frequency matched filter plus 
the following five filters in our study: 


The matched multi-filters (MMF) were first introduced by Her- 
ranz et al. (2002) for SZ detection. They were further developed 
by Melin et al. (2006) and extensively used to construct the suc¬ 
cessive Planck SZ cluster catalogs (Planck Collaboration et al. 
2011b, 2014a, 2015c). The MMF target the extraction of a single 
component across a set of maps, assuming a known frequency 
dependance and a spatial distribution. The filtered map returns a 
unbiased estimate, oa, of aA with minimal variance: 

dA = f d^x'V'/x) ■ m(x). (11) 


Similarly to the single frequency case, the MMF for the A emis¬ 
sion, 'Va (x), is expressed in Fourier space as 

'PXk) = cr2[a,,]p-'(k)-tXz,k). (12) 

The error on the amplitude is now given as 


cr[aq] = ^{dj} - {ciAf 


f 


d^k t*Uk)-P-'(k)-Q(z,k) 


- 1/2 


(13) 


with P(k) being the inter-band cross-power spectrum matrix with 
contributions from (non-d) sky signal and instrumental noise. It 
is the effective noise matrix for the MMF and, as for the single 
frequency power spectrum, can be estimated directly on the data, 
since the A signal (dust, synchrotron or tSZ) is small compared 
to other astrophysical signals over the sky patch. 

Introduced in Planck Collaboration et al. (2013) for a mix¬ 
ture of dust and tSZ signals, multi-component filtering deals 


- Single component MMF for dust; 

- Single component MMF for tSZ; 

- Two-component MMF for dust-ttSZ; 

- Two-component MMF for dust-tsynchrotron; 

- Three-component MMF for dust-HtSZ-i-synchrotron. 

3.3. From observed signals to physical QSO properties 

The amplitudes oa (determined with profiles cut at 0 < 505oo) are 
converted into spherically integrated quantities within R 500 by 

Aa - dA f dQ ta(0) X CaISRsoo ^ P 500 ), (18) 

*J9<595oq 

where the factor CaCSRsoo —» P 500 ) is the conversion from the 
volume within SRsoo to Rsoo, given the adopted profile. Similarly, 
we obtain the source flux density vector, F, as 

F = f r dQT^so(0)xC^so(5R5oo^P5oo) (19) 

•J0<505OO 

We express the quantities Adust, F, and Asynch in niJy and Atsz in 
arcmin^. 

Following Beelen et al. (2006) we estimate the total dust 
mass of the QSO from Adust as 

Mdust = Adust[857 GHz(l + z)] B 85 VGHz(i+z)[Fd], (20) 

with x(v) = xo ( 2498 GHz /''’ '^0 - 0.4g^'cm^ and Dl(z) the lumi¬ 
nosity distance. From the synchrotron flux density at 100 GHz, 
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we compute the monochromatic synchrotron luminosity accord¬ 
ing to 


r _ r 


( 21 ) 


C[<W)] = 


yc-'[(w)j 


-1 


(28) 


which is referenced to the rest-frame frequency of 100 GHz as¬ 
suming our power law with as. 

We define the intrinsic Compton parameter, Fsoo, as the 
quantity F 500 from Planck Collaboration et al. (2013), 

( 22 ) 

DfXz) being the angular distance to the QSO in Mpc. Applying 
Eq. (22) of Arnaud et al. (2010), we estimate the total mass of 
the halo as 


3.4.2. ^^-test for the hypothesis 

If the assumptions on the nature and number of components are 
correct, the measured flux densities, F,, and the reconstructed 
model flux densities, 2 (S i), (z), must be consistent within 
the uncertainties. We introduce the residual flux density for a 
given QSO as 

Ri^Fi-Y,A,{SA)i{z). (29) 

A 


P4 


/h\-' 

rA,szE(z)-"^'l 

\0.7j 

A,(z) 


1/ff 




(23) 


with a = 1.78 and Ax(z) = 2.925 x 10^^ x 0.6145 x (^) ‘ x 

(ifo^) arcmin^. The so-called mass bias parameter, 
b, accounts for any bias between the estimated mass and the true 
halo mass, Msoo.tme, for example from violation of hydrostatic 
equilibrium and/or X-ray instrument calibration. Planck SZ clus¬ 
ter counts require high values of the mass bias, ( 1 -h) ~ 0 . 6 - 0 .7, 
for consistency with the base ACDM model favored by Planck' 
measurements of the primary CMB anisotropies (see Planck 
Collaboration et al. 2014b, 2015e, and references therein). 

The M 500 - T 500 relation can be derived from the last two equa¬ 
tions; 


Msqo 

/ h 

T 500 

10'3 Mo 

\O. 7 j 

2.00 X 10“® arcmin^ 


CotTelation between the measured flux densities, F, and the 
model amplitudes. A, leads to the following expression for the 
variance of the residuals: 

C[R],v = {RiRj) - (RiXRj) 

= C[% - y (5,),. (z)C[A],,^ ( 5 ^), (z), (30) 

A,II 


where 

C[¥]ij = (FiFj) - {FiXFj) 



(f 

v qsoj 

.(k) 

tqso 

/k)P,y(k) 


n(k)P 

/./(k) 


(31) 


with dr 47rr^Tqso(r‘). The residuals Rj are expected 

to be Gaussian distributed with zero mean. We employ Eq. (27) 
and (28) with R, instead of w, to compute the average value, and 
define the for the average residual as 


3.4. Statistics 


;y2^(R)GC-'[<R)]-<R). 


(32) 


3.4.1. Average values for the full QSO sample 

As described before, the individual QSO signals are expected to 
be faint, so we need to average them over the population. Here, 
we detail the computation of this average. We adopt the generic 
notation vi>,i for A^, Mdust, Tsynch, T 500 or F,, arranging them into 
a vector w, and assume that they are Gaussian distributed. 

If these different components of w are not cotTelated, we em¬ 
ploy the usual inverse-variance weighted average: 


{wa) = 


Cr[<Wi)] = 


V 

1 

-i 

(ii'/l)p 

^ 0 
L p 

'-^[(WA)p\ 


^ crH(yVA)p\ 


V 

1 

- -1 


A 


o-2[(wd)p] 


(25) 


(26) 


The average residual vector (R) has seven components, one per 
frequency. We expect the;y^ to follow Student’s law with seven 
degrees of freedom, illustrated with simulations in Sect. 4.2. 
The value of depends on the number of components and 
their frequency dependance. In particular, the dust emission is 
parametrized by and ySj, so the x^ also depends on their val¬ 
ues. If we leave these two parameters free, the x^ would be ex¬ 
pected to follow Student’s law with 7-2-5 degrees-of-freedom. 

We study the average residual vector, (R), because it 
is the most useful quantity for identifying deviations be¬ 
tween observed signals and model predictions (dust, dust-HtSZ, 
dust-Hsynchrotron or dust-i-tSZ-i-synchrotron). Individual ob¬ 
tained from the residuals R, can also be computed, but given the 
large measurement uncertainties, study of the individual xj dis¬ 
tributions do not efficiently discriminate models. 

3.5. Marginaiization over dust properties 


with p the index over the QSO catalogue. 

In the correlated case, the average values are obtained ac¬ 
cording to 


<w) = 


yC-'[(w)^] 

. P 


Yj C ‘ [(w)p] ■ (w)p 


(27) 


The previous section considered computation of (vi>,i) and;y^ for 
hxed dust properties, described by the parameters and [3^. Al¬ 
though Adust is only weakly dependent on these parameters, the 
derived dust mass, Mdu.st (Eq. 20), is more sensitive to them. We 
use Powell’s algorithm to minimize the x^ relative to Td, [id to 
hnd their best-fit values. We also marginalize over Td, Pd when 
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determining the physical properties of the QSO environment, 

^dust5 ^synch ttnd 1^500' 

By “data” in the following, we will mean the combination of 
the Planck maps and the BOSS QSO positions. The (w^) being 
Gaussian distributed (Sec. 3.4.1), we write 


P(wA\Td,lidA3f'a) oc exp 


1 / wj - 

2 \ cr[<w,)] j 


(33) 


With f’(data|7’(j,j8d) exp^-y j, Bayes’ 


theorem gives 


P(rd,y6d|data) oc P{Ti,Pi) x exp 



(34) 


Assuming a flat prior on P(Td,j3d), we compute the probability 
of Wd, Td, and /5d given the data as, 

P(wa, 7’d,/3d|data) oc P(wA\Td,f^d, data) x B(7’d,y6d|data). (35) 


Averages are now calculated by marginalizing over Td and 
/3d and via an integration over vi>d: 


Wi 



waP(wa, Td,f3d\dsLtn)dTdd/3dWA, 


(36) 


with variance 


= JfJ 


(wa - WAfP(wA, Td,/3d\diitSi)dTdd/3ddwA- 


(37) 


4.2. Degrees of freedom 

We test the number of degrees-of-freedom for the;^'^ statistic de¬ 
fined in Sect. 3.4.2 by injecting into the Planck maps a mock 
catalogue with only dust emission from the QSOs and then ex¬ 
tracting the signal with the dust-only MMF at fixed Td - (7’d)input 
and ySd = (y6d)input- We bin our catalogue into 148 sub-catalogs of 
2000 QSOs each to compute 148 independent values. The 
cumulative distribution of these values is given in the top panel 
of Fig. 2 as the solid red line. Student’s cumulative distribution 
with varying degrees-of-freedom (dof) are shown as the dashed 
lines. The red line follows Student’s law with seven degrees-of- 
freedom, corresponding to the seven Planck maps. When leaving 
Td and jSd free and choosing their best-fit values as the point min¬ 
imizing the;^^^, we obtain the solid blue line with 7-2=5 degrees- 
of-freedom, as expected. 

We next inject a mock catalogue with both dust and tSZ 
emission and re-extract the QSO signals using a dust-HtSZ Al¬ 
ter to compute the;^f^. The results are shown in the bottom panel 
of Fig. 2. Fixing Td and /3d to the input values also leads to a 
X^ with seven dof. Leaving the two dust parameters free low¬ 
ers the x^ to five dof, as for the dust-only case. Although we 
are searching for an additional component in the data between 
the first and the second tests, the;^f^ conserves 7/5 dof for (Td, 
ySd) flxed/free respectively. This is because the correlations be¬ 
tween our residuals are taken into account in the covariance ma¬ 
trix C[(R)] in Eq. (32). In other words, C[(R)] changes when 
considering the dust-only MMF or the dust-HtSZ MMF, so the;^f^ 
statistic defined in Eq. (32) does not depend on the number of 
components assumed for the MMF 


4. Simulations 

We validate our methods by injecting simulated QSOs into the 
Planck maps and re-extracting their properties. The simulations 
are described in Sect. 4.1. Section 4.2 shows that our;^'^ statistic 
carries the expected number of degrees-of-freedom. We demon¬ 
strate that the MMFs properly recover the properties of the mock 
catalogs in Sect. 4.3, and in Sect. 4.4 we examine the sensitivity 
our results to the dust model. 

4.1. Mock catalogs and injected maps 

We build mock catalogs directly from the original BOSS sample, 
keeping the same QSO redshift distribution (Fig. 1) but draw¬ 
ing the sky positions at random outside the Planck point-source 
mask. We fix the QSO host halo mass M 500 = 2.7xlO'^Mo, close 
to the value found from the data in Sect. 5.4. We then compute 
the expected tSZ flux for each host from (z, M 500 ) using Eq. (23). 

We assign to each host halo a constant dust mass 
(4^du.st)input = 2.5 X 10^Mq emitting with a modified blackbody 
spectrum at (rd)input = 25 K, (J3d)mput - 2.5. When specified, we 
also implement variations in the dust temperature (7’d)input using 
a Gaussian distribution with mean (rd)input = 25 K and standard 
deviation (crx)i„put = 5K. 

With this model we compute the signal of each QSO in the 
Planck bands (from tSZ and dust) and we inject them directly 
into the maps using our source profile, Tqso, convolved with the 
individual channel beams. We artificially fix 6 ^ to 0.27 arcmin 
for the injected profile to avoid possible biases due to the mis¬ 
match between the extraction profile and the halo profile at low 
redshift. The sensitivity of our result to this assumption is stud¬ 
ied in Sect. 5.5. 


4.3. Extraction of simulated QSO 

We first consider null tests where the signals are extracted at ran¬ 
dom positions in the Planck maps. Table 1 shows the output for 
three Alters (dust, dust-ttSZ, dust-Hsynch) assuming fixed values 
of Td = 25K and /3d - 2.5 in the extraction. The output dust 
masses, SZ masses and synchrotron luminosities are compatible 
with zero, as expected, and the x^ do not show any significant 
deviation from Student’s law with the expected seven dof. 

Results of the extraction on mock catalogs injected into the 
Planck maps are summarized in Table 2. The recovered dust 
masses, SZ masses, and synchrotron luminosities are marginal¬ 
ized over Td and /3d as described in Sect. 3.5. 

Eor the injected dust-HtSZ catalogue, the dust-ttSZ Alter re¬ 
covers a unbiased estimate of the input dust temperature and 
spectral index. The input dust mass and the halo mass from the 
tSZ are also recovered at 5 /A ~ 8 and S /N ~ 14. The dust 
mass recovered using the dust-only Alter is biased by H-3.4cr due 
to contamination of the dust emission by the tSZ. The recov¬ 
ered dust spectral index is also biased by -3.7cr. Adopting the 
dust-Hsynch Alter leads to a signiflcant detection of negative syn¬ 
chrotron luminosity at H-4.6cr, which mimics the tSZ efl’ect at 
low frequencies. 

The dust-only extraction provides a reduced ;t'^/5 value of 
12.1, larger than the dust-HtSZ value of 1.1, indicating that the 
dust-HtSZ model must be preferred over the dust-only model. The 
dust-Hsynch model provides a better;^f^/5 than the dust-only case 
(4.0), but the signiflcantly negative value for the synchrotron lu¬ 
minosity discards the model. The (Td,/3d) contours for the three 
Alters are displayed in Fig. 3. The dust-ttSZ Alter shows contours 
enclosing the input values while the two other Alters are biased. 

We also implemented the triple dust-HtSZ-i-synch Alter and 
applied it to the dust-HtSZ simulation. The triple Alter trans- 
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Fig. 2. Cumulative distributions when I’d and Pi are fixed in the fil¬ 
tering (solid red curve) or when they are left free (solid blue curve). Top 
panel: Distribution for a dust-only filtering on a dust-only mock cata¬ 
logue. Bottom panel: Distribution for a dust-l-tSZ filtering on a dust-l-tSZ 
mock catalogue. In both cases (dust-only or dust-HSZ), fixing Ti and pi 
leads to a x^ distribution with seven degrees-of-freedom (correspond¬ 
ing to our seven maps), while leaving them free reduces the number of 
degrees-of-freedom to five. 


forms part of the tSZ emission into a negative (i.e., unphysical) 
synchrotron emission, as shown in Table 2. It recovers biased 
Tj, Pi, and The reduced x^ 1^ increases with respect to 

the dustH-tSZ filter, showing that the dust-HtSZ fit must be pre¬ 
ferred over the dust-HtSZ-Hsynch fit. We show in Sect. 5 that the 
triple filter exhibits the same behavior on the data. We will thus 
present our main results using the dust, the dustH-tSZ, and the 
dust-Hsynch filters, depending on the value of the reduced 

4.4. Possible systematics due to the dust characteristics 

We test the sensitivity of the extraction method to improper mod¬ 
eling of the dust by injecting mock QSOs with a dust tempera¬ 
ture following a Gaussian distribution, described in Sect. 4.1 , and 
then recovering the parameters with the different MMF. Results 
are shown in the bottom part of Table 2. The dust (ctj = vK)-i-tSZ 
filters are a dust-HtSZ filter for which a Gaussian scatter of x 
Kelvin in the dust temperature is added. As expected, the dust 
(cTj -5 K)-i-tSZ filter provides an unbiased estimate of the pa¬ 
rameters. The derived dust and SZ masses are not very sensitive 


Fig. 3. Contours of x^ for the mock QSO dust+tSZ catalogue using 
the dust, the dustr-tSZ, and the dust-l-synch filters. The solid black line 
indicates the Icr deviation from the minimal x^ and the dotted black 
line the 2cr deviation. The horizontal and vertical dashed lines mark the 
position of the input Td and pi. 


Of = 2 K CTt = 5K aT=10K 



2 3 2 3 2 3 

Pi Pi Pi 


Fig. 4. Contours ofx^ for the mock QSO dust (ctt^S K)-l-tSZ catalogue, 
obtained using the dustfiTT =2 K)-l-tSZ, the dust (ctj =5 K)-l-tSZ, and the 
dust (cr-p =10K)-l-tSZ filters. The solid black line indicates the Icr devi¬ 
ation from the minimal and the dotted black line the 2cr deviation. 
The horizontal and vertical dashed lines mark the position of the input 
Td andjSd. 


to the dispersion in the dust temperature; the dustH-tSZ filter re¬ 
turns a satisfactory estimate of the two quantities. However, for 
this latter filter the recovered dust parameters and Pi are sig¬ 
nificantly biased. This sensitivity is illustrated in Fig. 4 display¬ 
ing the iTi,Pi) contours for the three last filters of Table 2. 


5. Results 

We apply the method described in Sect. 3 to the real Planck data. 
In Sect. 5.1, we show the average total flux of the QSO sample 
across the Planck channels. We then estimate the contribution 
of the various components in Sect. 5.2 (dust), 5.3 (synchrotron), 
and 5.4 (hot gas). 
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Table 1. Null test performed at fixed (7^=25 K, p^=2.5). The filtering is done at random positions in the Planck maps. 


Filter 

x^n 

Ti 

(K) 

Pa 

(^dust) 

(10'Mo) 

<F500> 

(10“® arcmin^) 

{^synch) 

Mock catalogue : Null test 
input 

_ 

0 

0 

0 

dust 

0.7 

25 

2.5 

0.003 + 0.012 

- 

- 

dustH-tSZ 

1.2 

25 

2.5 

0.003 + 0.012 

0.77 + 0.78 

- 

dust-i-synch 

3.5 

25 

2.5 

0.003 + 0.012 

- 

-0.05 + 0.13 


Table 2. Extracted QSO properties for the mock catalogs, averaged over (Td./ld) as described in Sec. 3.5. The first mock catalogue (upper part of 
the Table) features both the tSZ and the dust emission. The second catalogue (bottom part) also includes variation in the dust temperature Tj. For 
each case the input line gives the injected values. The other lines shows the values extracted with the various filters. For each simulation, the filter 
assuming the same components as the input is highlighted in gray. 


Filter 

X^/5 

Ti 

(K) 

Pa 

^dust 
(10'Mo) 

Z500 

(10“® arcmin^) 

M 500 

(10''Mo) 

^synch 

Mock catalogue ; dust+tSZ 
input 

25 

2.5 

2.50 

12.02 

2.74 

0 

dust 

12.1 

26.1+0.6 

2.24 + 0.07 

3.54 + 0.31 

- 

- 

- 

dust+tSZ 

1.1 

25.4 + 0.8 

2.43 + 0.1 

2.71+0.34 

13.03 + 1.63 

2.86 + 0.20 

- 

dust+synch 

4.0 

27.4 + 0.8 

1.95 + 0.08 

5.18 + 0.42 

- 

- 

-0.65 + 0.14 

dust+tSZ+synch 

2.9 

27.4+ 1.0 

2.05 + 0.13 

4.29 + 0.65 

8.68+1.90 

2.27 + 0.30 

-0.37 + 0.17 

Mock catalogue ; dust+tSZ with Ti Gaussian distributed 
input - 25 + 5 2.5 

2.50 

12.02 

2.74 

0 

dust 

12.7 

30.9+1.2 

2.06 + 0.08 

3.27 + 0.22 

- 

- 

- 

dust+tSZ 

1.1 

29.8 + 0.8 

2.27 + 0.08 

2.46 + 0.27 

13.03 + 1.70 

2.86 + 0.21 

- 

dust+synch 

3.7 

33.6 + 0.4 

1.75 + 0.05 

4.69 + 0.29 

- 

- 

-0.69 + 0.13 

dust(cr7-=2 K)+tSZ 

1.1 

29.3 + 0.8 

2.28 + 0.08 

2.47 + 0.27 

12.81 + 1.68 

2.83 + 0.21 

- 

dust(cr7-=5 K)+tSZ 

1.1 

26.4+ 1.2 

2.37 + 0.10 

2.67 + 0.29 

12.97+1.68 

2.85 + 0.21 

- 

dust(o-7-=10 K)+tSZ 

1.2 

7.5 + 1.6 

2.95 + 0.08 

5.43 + 0.40 

11.45 + 1.31 

2.66 + 0.17 

- 


5.1. The QSO signal at Planck frequencies 

We first estimate the mean flux density, (F), of the 291 165 QSOs 
selected in Sect. 2.1 by averaging the individual flux densities in 
each band as extracted with a single frequency matched Alter 
centered on the QSO positions. The average total flux is shown 
as the black diamonds in Fig. 5 (main panel and inset). It is 
strictly positive in all Planck bands above 100 GHz, and positive 
but compatible with zero at 70 GHz. The signal continuously in¬ 
creases from the lowest to the highest frequencies; emission in 
the direction of the QSO is dominated by dust. Since none of the 
frequencies below 217 GHz presents a negative flux, there is no 
obvious indication of strong tSZ emission. 

Using the dust-only MMF with and fii adjusted to min¬ 
imize the;^^^ defined in Eq. (32) (Td = 18.7 K, Pi - 2.79), we 
compute the mean dust flux density Adu.st and plot the resulting 
mean dust spectrum as the red triangles in Fig. 5. Dust accounts 
for essentially all the observed signal. The residual (difference 
between the black diamonds and the red triangles) is shown in 
the bottom panel of the figure, and is consistent with zero ex¬ 
cept at 353 GHz. This deviation explains why the reduced 1^ 
highlighted in the first gray line of Table 3 is not good ix^ jS- 
5.2). 

By marginalizing over and p^, we obtain an average dust 
mass, Mdust = (0.84 + 0.07) x IO^Mq. This value is ~ 25 times 
smaller than the estimated dust mass in optically-selected clus¬ 
ters, around 2 x 10 ®Mq (Gutierrez &l Lopez-Corredoira 2014). 
It is also five times smaller than the lowest estimated QSO dust 
mass of the Beelen et al. (2006) sample. We note, however, that 
our average mass is dominated by low redshift objects (z ~ 0.5), 
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Fig. 5. Average flux density for the QSO population across Planck chan¬ 
nels (black diamonds), and dust signal estimated with the dust-only fil¬ 
ter (red triangles) for and at the minimum The bottom panel 
shows the average residuals. 
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as shown in Fig. 7, while the QSO sample studied in Beelen et al. 
(2006) resides at higher redshift. 

The recovered dust mass also depends strongly on the as¬ 
sumed value for and p^. If we fix Pi — 1.6, as in Beelen et al. 
(2006), we find Mdust = (1.70 + 0.08) x IO^Mq, in better agree¬ 
ment with the above-cited study, but the;t'^/5 increases to 16.3, 
as shown in Table 3. This result demonstrates that the bulk of 
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the BOSS QSO population contains significantly less dust than 
QSOs examined in these previous studies. 

For the dust temperature, we find a marginalized — 19.1 + 
0.8 K, close to the dust temperature of normal galaxies (e.g., 
Clemens et al. 2013). This is significantly lower than the tem¬ 
perature found by Beelen et al. (2006) and Dai et al. (2012), 
ranging between 33 and 55 K, and 18.1 and 79.7 K respectively. 
We note that Dai et al. (2012) use the two component model 
from Blain et al. (2003) with Pd — 2, and Beelen et al. (2006) fix 
Pd — 1.6 for the majority of their QSOs. We find ySd = 2.71+0.13, 
significantly higher than these values. When fixing Pd — 1.6 in 
our analysis, our result moves along the degeneracy line in the 
(T(i,j0d)-plane to reach Td - 28.2 + 0.5 K, in better agreement 
with the values published in the previously mentioned analyses. 
Restricting this test to high redshift QSOs leaves our conclusions 
unchanged (line 6 of Table 3). 


5.2. Redshift dependence of the dust properties 

In order to adapt the MMF to a possible variation of the QSO 
dust properties across redshift, we divide our QSO sample into 
nine regular bins of size Az = 0.5, shown as the vertical dashed 
lines in Fig. 1. We include in the last bin the QSOs with 4 < z < 5 
because there are not enough statistics to fill a tenth bin. 

We minimize the for each bin individually and plot the 
value for the combination Td x pf^^ in Fig. 6. This quantity fol¬ 
lows the degeneracy line for the dust parameters. There is no 
evidence for significant evolution across redshift between z = 0 
and 2. For z > 2, Td x pf^^ increases with redshift. If we suppose 
that Pd remains constant, this implies that the dust was hotter 
at high redshift and cooled until z = 2 before stabilizing at the 
current values. 

The dust flux density extracted using the dust-only MMF 
is shown in Fig. 7 as the black diamonds. It varies with red¬ 
shift, increasing from z = 0toz = 2(~8 mJy within R 500 at 
857 GHz) and decreasing with redshift for z > 2. The corre¬ 
sponding dust mass (Eq. 20) is shown by the blue diamonds and 
follows the same trend as the flux density, reaching a maximum 
of ~ 5 X IO^Mq at z ~ 2. This is of the same order as dust masses 
determined from high signal-to-noise observations of individual 
QSOs (e.g., Beelen et al. 2006; Wang et al. 2007). 

The IR luminosity, or equivalently the dust mass (Eq. 20), is 
a tracer of star formation, and it is remarkable that the dust mass 
evolution in Eig. 7 follows that of the general star formation rate 
with a peak around z ~ 2 (See, e.g., Eig. 9 in Madau & Dickinson 
2014). The accretion rate onto AGN SMBHs follows the same 
trend (Hopkins et al. 2007; Shankar et al. 2009; Aird et al. 2010; 
Delvecchio et al. 2014). This is strong evidence that star forma¬ 
tion in QSO environments is typical of the general galaxy popu¬ 
lation, and also that it is linked to the QSO central engines. Our 
result clearly shows this trend for a large, representative sample 
of QSOs, and is consistent with the study by Wang et al. (2015) 
of the correlation between QSOs and Herschel measurements of 
the cosmic infrared background. 

Table 3 also contains results for the dust+tSZ, dust+synch, 
and the dust+tSZ+synch filters on the full QSO population 
(0.1 < z < 6.5). All three of the multiple component filters in¬ 
crease the with respect to the dust-only filter. 
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Fig. 6. Top panel: Average T;j as a function of redshift z. Middle panel: 
Average Pd as a function of redshift z. There is a clear evolution between 
the low-z (z<1.5) , high-temperature and low-y6 and the high-z, low- 
temperature and high-j8 QSO populations. Bottom panel: Average dust 
Td X /f} ® (combination describing the degeneracy line) as a function 
of redshift z. There is no evidence for evolution of the dust properties 
between z = 0 and 1.5. The product increases between z = 1.5 and 4. 
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Table 3. Average properties of different QSO samples in Planck data, marginalized over and (see Sec. 3.5), except for line 1 and 6 marginal¬ 
ized over Td only, being fixed to 1.6 as in Beelen et al. (2006). QSO with at least one radio counterpart at 1.4GHz are flagged with ’yes’ in the 
FIRST column. QSO with no identified counterpart are flagged with ’no’. For each sample, the filter leading to the smallest is highlighted in 
gray. The different sample sizes are as follows: full sample (0.1 < z < 5), 291 165; restricted redshift (2.5 < z < 4), 90 576 of which 3354 have 
FIRST counterparts and 77 665 do not, and the remaining 9557 fall outside the FIRST coverage. 


Zmin 

Zmax 

FIRST 

Filter 

X^I5 

Ti 

(K) 

Pi 

^dust 
(lo** Mo) 

5^500 

(10“^ arcmin^) 

Msoo 

(10*3 Mo) 

^synch 

(10-3 

0.1 

5 

- 

dust 

16.3 

28.2 ±0.5 

1.6 (fixed) 

1.70 ±0.08 

- 

- 

- 

0.1 

5 

- 

dust 

5.2 

19.1 ±0.8 

2.71 ±0.13 

0.84 ± 0.07 

- 

- 

- 

0.1 

5 

- 

dust-rtSZ 

9.1 

18.6 ±0.8 

2.80 ±0.13 

0.76 ± 0.07 

17.02 ±0.82 

3.32 ± 0.09 

- 

0.1 

5 

- 

dust+synch 

6.7 

20.0 ± 1.0 

2.52 ±0.14 

0.98 ± 0.09 

- 

- 

0.34 ±0.13 

0.1 

5 

- 

dust+tSZ+synch 

16.1 

21.3± 1.1 

2.24 ±0.15 

1.21 ±0.12 

17.15 ±0.82 

3.34 ± 0.09 

0.66 ±0.13 

2.5 

4 

- 

dust 

12.4 

28.7 ±0.5 

1.6 (fixed) 

7.85 ± 0.42 

- 

- 

- 

2.5 

4 

- 

dust 

4.2 

21.9 ±0.8 

2.63 ±0.14 

2.82 ± 0.46 

- 

- 

- 

2.5 

4 

- 

dust+tSZ 

2.3 

22.5 ± 0.9 

2.59 ±0.18 

2.60 ± 0.53 

8.62 ± 1.36 

2.26 ± 0.20 

- 

2.5 

4 

- 

dust+synch 

3.8 

23.8 ± 0.9 

2.25 ±0.16 

4.46 ± 0.82 

- 

- 

-26.41 ±5.78 

2.5 

4 

- 

dust+tSZ+synch 

3.4 

23.5 ± 1.0 

2.34 ± 0.22 

3.75 ± 1.05 

5.92 ± 2.25 

1.63 ±0.71 

-11.37 ±8.66 

2.5 

4 

yes 

dust 

4.9 

32.5 ±5.1 

1.24 ±0.44 

16.31 ±6.40 

- 

- 

- 

2.5 

4 

yes 

dust-HtSZ 

2.0 

23.3 ±3.2 

2.01 ± 0.46 

14.32 ± 6.24 

-40.14 ±7.57 

- 

- 

2.5 

4 

yes 

dust+synch 

0.4 

18.8 ±2.8 

3.72 ± 0.67 

1.11 ± 1.14 

- 

- 

236.91 ± 30.47 

2.5 

4 

yes 

dust+tSZ+synch 

0.3 

15.7 ±3.3 

5.33 ± 1.43 

0.38 ± 0.68 

15.77 ± 10.44 

2.23 ± 1.88 

278.76 ± 40.69 

2.5 

4 

no 

dust 

4.4 

21.8 ±0.9 

2.69 ±0.19 

2.57 ± 0.54 

- 

- 

- 

2.5 

4 

no 

dust+tSZ 

1.5 

22.4 ± 1.2 

2.68 ± 0.23 

2.18 ±0.56 

10.83 ± 1.46 

2.58 ±0.20 

- 

2.5 

4 

no 

dust+synch 

2.7 

24.0 ± 1.0 

2.22 ±0.18 

4.47 ± 0.89 

- 

- 

-34.17 ±6.23 

2.5 

4 

no 

dust+tSZ+synch 

2.3 

23.6 ± 1.1 

2.34 ± 0.23 

3.65 ± 1.09 

6.82 ±2.41 

1.83 ±0.67 

-16.86 ±9.24 

2.5 

4 

no 

dust(crx=2 K)+tSZ 

2.1 

20.9 ± 1.1 

2.81 ± 0.22 

2.03 ± 0.45 

10.55 ± 1.43 

2.54 ± 0.20 

- 

2.5 

4 

no 

dust(crT=5 K)-HSZ 

1.9 

10.7 ± 1.5 

3.61 ± 0.23 

3.08 ± 0.33 

9.37 ± 1.49 

2.37 ±0.21 

- 

2.5 

4 

no 

dust(crT = 10 K)-HSZ 

11.7 

6.6 ±0.9 

2.31 ±0.02 

13.07 ± 1.12 

9.73 ± 1.47 

2.42 ± 0.21 

- 



z 


Fig. 7. Average dust flux density after marginalization over T^ and pi 
(black diamonds) as a function of redshift z. The corresponding dust 
mass is shown as blue triangles (Eq. 20). The extraction is achieved 
using the dust-only filter. 


5.3. Synchrotron emission 

We now select only QSOs with at least one FIRST^ radio source 
counterpart at 1.4 GHz using the FIRST_MATCH keyword in 
the DR12 catalogue (FIRST_MATCH=1). This provides a sub¬ 
catalogue of 9983 QSOs. Applying the dust-Hsynch filter to this 
sub-catalogue using the same redshift binning, we plot the av¬ 
erage dust and synchrotron flux densities, Adust and Asynch, in 
Fig. 8. Dust emission from the sub-catalogue is compatible with 
that of the full catalogue, although with larger uncertainties due 
to the smaller sample (blue diamonds). 

This population of radio-loud QSOs exhibits strong syn¬ 
chrotron emission, as shown in the lower panel. The monochro- 


^ http://sundog.stsci.edu 


matic synchrotron luminosity (green diamonds) increases from 
0 to reach 0.2 Lq Hz * at the higher redshifts; We note that this 
is the luminosity at the rest-frame frequency 100 GHz, assuming 
our power-law in Eq. (eqilsynch). 

Fig. 9 compares the synchrotron emission determined from 
Planck (expressed at 100 GHz) to the FIRST measurements at 
1.4 GHz. For this comparison, we bin-average the FIRST sig¬ 
nal with the same MMF weights used to extract the synchrotron 
signal. The emission seen by Planck steadily decreases with red¬ 
shift relative to the signal measured by FIRST, demonstrating a 
variation in synchrotron index with redshift. 

We show the effective spectral index (for an assumed power 
law) between FIRST and Planck in Fig. 10. The spectrum is 
quite flat at low redshift with values for the spectral index that 
are much smaller than our adopted value of 0.7. However, these 
are values describing the emission over a very large frequency 
range, while our adopted value is assumed to describe the emis¬ 
sion around 100 GHz. 

We also see that the effective spectral index steepens with 
redshift. This could be due to evolution intrinsic to the sources, 
or more simply the effect of curvature in the synchrotron spec¬ 
trum. It is important to note in this light that although our syn¬ 
chrotron luminosity measurements (Eq. 21) are referenced to 
100 GHz rest-frame, assuming our adopted power law, the ef¬ 
fective index shown in Eig. 10 is taken between the observed 1.4 
and 100 GHz bands. In other words, it is the effective index be¬ 
tween rest-frame frequencies of 1.4(1 + z) and 100(1 -H z)GHz. 
The apparent evolution with redshift in the figure could therefore 
be steepening in the synchrotron emission with frequency that 
would be expected given the greater energy loss of the higher 
energy electrons contributing to the signal as we move up in rest- 
frame frequency. 

5.4. Hot halo gas 

The strong synchrotron emission of QSOs with EIRST counter¬ 
parts complicates extraction of the tSZ signal from the hot gas 
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Fig. 8. Top panel: Average dust flux density (black diamonds) for QSOs 
with a radio counterpart in the FIRST catalogue. The signal is extracted 
with the dust+synch filter and marginalized over Tj and ySj. The corre¬ 
sponding dust mass is shown as the blue triangles. Bottom panel: Av¬ 
erage synchrotron flux density at 100 GHz observed frequency (black 
diamonds) and the corresponding monochromatic luminosity (magenta 
triangles). 




U 
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Fig. 9. Average FIRST flux density at 1.4 GHz observed frequency com¬ 
pared to the average synchrotron flux density at 100 GHz observed fre¬ 
quency as a function of redshift z. The magenta triangles correspond to 
the black diamonds in the bottom panel of Fig. 8. 


in the host halo. We therefore construct a sub-sample without 
FIRST counterparts (FIRST_MATCH=0) of 252 888 QSOs and 
apply the dust-HtSZ filter. Results are shown in Fig. 11, and sum¬ 
marized in the bottom segment of Tab. 3 for this and other filters. 

Applied to the same redshift range, the other filters 
yield larger values. In particular, the dust-i-synch and 
dust-i-tSZ-i-synch filters show the same behavior as on simu¬ 
lations; they convert part of the tSZ signal into negative syn¬ 
chrotron luminosity. The dust((T 7 - =x K)-i-tSZ filters presented in 
Sect. 4.4 also increase the favoring a low scatter in the dust 
temperature. 

The filter does not detect the tSZ effect at z < 2.5. The tSZ 
signal at z < 1.5 is slightly negative at ~ 2 cr, pointing towards a 
possible contamination of the extracted signal at low z by a weak 
synchrotron emission. The filter does show a clear signal rising 
with redshift over the range 2.5 < z < 4. The global significance 


Fig. 10. Effective spectral index of power law in frequency between 
1.4 GHz (FIRST) and 100 GHz {Planck) observation frequencies. 


of the signal in the latter redshift range is 7.4cr when expressed 
in terms of the population mean, Tsoo. 

The dust-HtSZ-Hsynch filter applied on the radio-loud sub¬ 
catalogue, i.e., with at least one FIRST source (third segment 
of Tab. 3), provides a slightly better fit (t^/ 5 = 0.3) than the 
dust-Hsynch filter (t^/ 5 = 0.4) over the 2.5 < z < 4 range. 
The recovered SZ signal remains compatible with the value de¬ 
rived from the radio-quiet sub-catalogue, but the uncertainties 
are large with the signal appearing at only 1.5cr. 

When applying the tSZ-only filter to the 2.5<z<4 popula¬ 
tion, we find SZ-based masses that are twice as large as the mass 
provided by the dustH-tSZ filter; the fit, however, is quite bad at 
n > 100 (We note that there are seven degrees-of-freedom 
in this case because the dust temperature and emissivity do not 
come into play). This simply reflects the fact that a tSZ-only sig¬ 
nal is a poor description of the spectrum of the source. Details of 
these results are given in Table 4. 

We stacked the ROSAT All-Sky Survey maps'^ of the radio¬ 
quiet QSO sub-sample to find a significant signal, correspond¬ 
ing to a population mean luminosity of L 500 = (3.30 + 0.30) x 
10'^'*erg/s in the [0.1-2.4] keV band, the error taken as the stan¬ 
dard deviation of the stacked maps. Pratt et al. (2009) estab¬ 
lished a power-law scaling relation between X-ray luminos¬ 
ity and halo mass based on analysis of a representative sam¬ 
ple of X-ray galaxy clusters employing Bivariate Correlated Er¬ 
rors and intrinsic Scatter estimators with a minimization of the 
residuals in X-ray luminosity (BCES (YIX)). If we adopt their 
Malmquist-bias corrected relation between L 500 , in this energy 
band, and mass (line 6 of their Table B.2), we infer a mean 
halo mass of M 500 = 5.87 x IO'^Mq. This is 2.3 times higher 
than our SZ-based mass; alternatively, it implies an expected 
T 500 = (51.4 + 4.6) X lO^^arcmin^, or 4.7 times larger than our 
measurement. In other words, the tSZ and X-ray signals do not 
follow the empirical relation for thermal halo gas emission seen 
at low redshift (Planck Collaboration et al. 2011c). This is not 
a surprise given that we expect a non-negligible contribution to 
the X-ray emission from the central engine of the QSO itself. 


http: //cade. irap. omp. eu/dokuwiki/doku. php?id=rass 
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Table 4. Average properties from a SZ filtering for several samples of QSO. 


Zmin 

Zmax 

EIRST 

Eilter 

x^n 

(Tsoo) 

(10^^ arcmin^) 

0.1 

5 

- 

tSZ 

520.8 

20.25 + 0.81 

2.5 

4 

- 

tsz 

141.2 

29.94+ 1.17 

2.5 

4 

yes 

tsz 

17.1 

7.46 + 5.84 

2.5 

4 

no 

tsz 

119.7 

30.53 + 1.25 



z 


Fig. 11. Top panel: Average dust flux density (black diamonds) of the 
radio-quite QSO subsample, i.e., without a radio counterpart in the 
FIRST catalogue. The signal is obtained with the dust+tSZ filter and 
is marginalized over Tj and Pi. The corresponding dust mass is shown 
by the blue triangles. Bottom panel: Average Compton parameter (black 
diamonds) and corresponding integrated Compton parameter (red trian¬ 
gles) as a function of redshift. 

5.5. Sensitivity to the assumed size fls 

We fix fls = 0.27 arcmin throughout our analysis. The value cor¬ 
responds to a halo mass of M 500 = lO'^ Mq at z = 2. To test the 
robustness of our results to this assumption, we fix alternatively 
the size to 0s=O. 12 and 0s=O.57 arcmin, corresponding to halos 
of mass M 500 = 10 '^ Mq and M 500 = IO'^^Mq, respectively, at 
z-2. Results are shown in Tab. 5. Changing the filtering scale 
slightly increases the and only weakly impacts the derived 
parameters (< O.Scr), leaving our results essentially unchanged. 


6. Discussion 

Submillimeter observations are a valuable source of information 
on the physical conditions in QSO environments because they 
probe star formation in and around the hosts and energy injec¬ 
tion into the surrounding medium by the SMBHs. Due to atmo¬ 
spheric absorption, this waveband is however difficult to access 
from the ground, and this has limited studies to relatively small 
samples of several tens of objects (Omont et al. 2001; Isaak et al. 
2002; Omont et al. 2013; Priddey et al. 2003; Beelen et al. 2006; 
Wang et al. 2007; Omont et al. 2013). 

The space-based platforms WMAP, Planck and Herschel 
now open this window for much more extensive studies. Tak¬ 
ing advantage of Planck's all-sky and wide-frequency coverage, 
we have determined several mean characterstics of the QSO pop¬ 
ulation through binning measurements of dust, synchrotron and 
tSZ signals from the very large BOSS QSO sample. 


6.1. Dust emission 

We find that the QSO spectrum is dominated by graybody dust 
emission with a temperature of - 19.1 + 0.8 K and an emis- 
sivity spectral index of = 2.71 +0.13 when averaged over the 
entire sample spanning the full redshift range, 0.1 < z < 5; these 
are the one-dimensional marginalized constraints for the dust- 
only filter that gives the best x^ (see Tab. 3). While there is no 
clear evidence that these quantities evolve with redshift at z < 2, 
there is a rise at higher redshifts in the product describing 

the degeneracy direction, at higher redshifts (Fig. 6). 

Curiously, the temperature is notably lower than found in 
ground-based observations of high-redshift quasars, which find 
temperatures of 40-50 K, while our value is more typical of or¬ 
dinary galaxies like the Milk Way. Our dust emissivity index is 
significantly steeper than the values of 1-2 for normal Galac¬ 
tic dust, and should be compared to the value of 1.6 found in 
the ground-based surveys. It is notably steeper than the limiting 
case of p - 2 for ideal insulators and conductors far from a reso¬ 
nance. These characteristics remain in the other sub-samples and 
component filters summarized in Tab. 3, apart the unusual be¬ 
haviour of the radio-loud sub-sample with FIRST counterparts. 
Fixing Pd — \ .6, we obtain a somewhat higher mean tempera¬ 
ture, Tdust = 28.2+0.5 K, but at the cost of a much worse spectral 
fit according to the large increase inx^ (see Tab. 3). 

We calculate dust mass using the simple prescription of 
Eq. (20), which can also be viewed as a measure of far-IR lu¬ 
minosity. An important result is our finding that the dust mass, 
a tracer of local star formation, evolves with redshift in the 
same way as the global cosmic star formation rate, peaking at 
z ~ 2 (See e.g. Fig. 9 in Madau & Dickinson 2014). Star forma¬ 
tion in QSO environments apparently does not distinguish itself 
from others. It is also known, from measures of AGN luminosity 
functions, that the SMBH accretion rate follows the same evo¬ 
lutionary track (Hopkins et al. 2007; Shankar et al. 2009; Aird 
et al. 2010; Delvecchio et al. 2014). 


6.2. Synchrotron emission 

A subsample of our QSOs are radio loud with counterparts at 
1.4 GHz in the FIRST catalogue. We detect synchrotron emis¬ 
sion from this subsample with an monochromatic luminos¬ 
ity, referenced to rest frame frequency 100 GHz (adopting our 
power-law of Eq. 21), increasing with redshift between z = 0 
and z = 3. Comparing the EIRST and Planck flux densities, we 
find an effective spectral index between the two corresponding 
rest-frame frequencies that increases with redshift from as = 0.3 
to as = 0.7. This steepening is likely the result of curvature in 
the spectrum as we observe progressively higher rest-frame fre¬ 
quencies. 
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Table 5. Average properties of the 2.5 < z < 4 QSO sample without FIRST counterpart for three filtering sizes 0^. 


Zmin 

Zmax 

EIRST 

0 

(arc min) 

Eilter 

a"/5 

Td 

(K) 

Pi 

^dust 

(10*^ Mo) 

Z500 

(10“® arcmin^) 

M500 

(lO'^Mo) 

2.5 

4 

no 

0.12 

dust-HtSZ 

2.0 

22.3+ 1.0 

2.68 + 0.21 

2.00 + 0.47 

10.89+ 1.42 

2.59 + 0.19 

2.5 

4 

no 

0.27 

dust-HtSZ 

1.5 

22.4+ 1.2 

2.68 + 0.23 

2.18 + 0.56 

10.83 + 1.46 

2.58 + 0.20 

2.5 

4 

no 

0.57 

dust-HtSZ 

2.6 

22.6+ 1.0 

2.67 + 0.20 

2.23 + 0.50 

11.15 + 1.47 

2.62 + 0.20 


6.3. tSZ and feedback 

Restricting analysis to the radio-quite sub-sample, we detect the 
tSZ effect (with the dust-rtSZ filter) at 5/A^ ~ 7.4 over the red- 
shift range 2.5 < z < 4, with an integrated Compton-y param¬ 
eter of F 500 = (10.83 + 1.46) X 10“®arcmin^ (one-dimensional, 
marginalized over dust properties). We do not, however, find ev¬ 
idence of tSZ signal at the lower redshifts. At z < 1.5, we obtain 
a mean value of ? 5 oo(z < 1.5) = (-4.93 + 2.48) x 10“® arcmin^, 
which we convert simply into an upper limit at two sigma above 
zero as 

?500(z < 1.5) < 4.97 X 10^® arcmin^ (2o-). (38) 

Similarly, we find an upper limit at z < 2.5 of 

FsooCz < 2.5) < 2.93 x 10^® arcmin^ (2o-), (39) 

even stronger, as expected given the behaviour seen in Fig. 11. 

Ruan et al. (2015) recently reported a detection of the tSZ 
signal by stacking the Hill & Spergel (2014) tSZ map, con¬ 
structed from the Planck frequency maps, on 26 686 spectro¬ 
scopic DR7 QSOs. For their low redshift bin at z < 1.5, they 
find Fjqq = (74 ± 30) x 10“® arcmin^, quoted in terms of the 
Compton-y parameter integrated along the line-of-sight cylin¬ 
der. The conversion factor between the cylindrical and spheri¬ 
cal integrated Compton-y parameters for our adopted profile is 
Fjoo - 1-52 X ? 5 oo, which translates the upper limit of Eq. (38) 
to < 7.55 X 10“® arcmin^ (2cr), well below their detection 
level. 

One possible explanation for the difference is dust contam¬ 
ination. Ruan et al. (2015) took particular pains to correct for 
any contamination by dust in the tSZ map. Our filters directly 
measure the individual local signals from dust and tSZ, as well 
as any covariance between them, allowing us to separate and 
marginalize over the dust influence when quoting our tSZ sig¬ 
nal strengths. The importance of separating out the dust signal is 
easily appreciated from the results of applying the tSZ-only filter 
given in Tab. 4. Contamination by dust in the tSZ-only filter can 
produce tSZ signals ~ 3 times larger. 

The tSZ signal directly measures the thermal energy stored in 
diffuse ionized gas in the QSO environment. We expect a signal 
from the intra-halo gas at the virial temperature of the QSO host 
halo, and there may be additional heating powered by feedback 
from the SMBH. Chatteijee et al. (2010) and Ruan et al. (2015) 
both interpreted their tSZ signals as evidence for feedback into 
the surrounding gas. 

To look for this, we estimate the tSZ signal expected from 
intra-halo gas using the scaling relation between halo mass and 
tSZ signal in Eq. (23). We first establish a mean mass scale for 
the sample of 77 655 QSOs in 2.5 < z < 4 over the 14 555 sq. 
deg. SDSS footprint. Using the Tinker et al. (2008) mass func¬ 
tion, we find the lower mass bound above which there are as 
many predicted halos as observed QSOs in this redshift range 
and for this sky area. We then calculate the mean mass by inte¬ 
grating over the mass function above this mass limit and over this 


redshift range. This gives us an estimate of the maximum mean 
mass characterizing the sample; QSOs could be hosted by lower 
mass halos if the QSO selection function is not 100% complete, 
as is most certainly the case. 

We find a mean halo mass of Msoo.tme = 2.18 x 10'^ Mq, 
a value in reasonable agreement with Richardson et al. (2012), 
who found a QSO halo mass 14.Ugg x IO^^/i^'Mq at z ~ 3.2 
for the SDSS sample. Our characteristic mass corresponds to a 
predicted tSZ signal of T 500 = 8.47(1 - f?)'-’* x lO^^arcmin^, 
only one sigma below our measured signal of >500 = (10.02 + 
1.34) X 10“® arcmin^ if b - 0. This would not leave much room 
for additional heating by feedback. With the higher mass bias 
values suggested by Planck" cluster counts, e.g., - b) - 0.6, 

we instead have Tsoo = 3.4 x 10“^ arcmin^, implying that ~ 60% 
of the gas’s thermal energy is generated by feedback. 

X-ray emission from the hot gas is also expected. Stack¬ 
ing the RASS maps at the QSO position, we find a luminosity 
L 500 = (3.30 ± 0.30) X lO^^erg/s in the [0.1-2.4]keV band. The 
tSZ - X-ray luminosity scaling relation established by Planck 
Collaboration et al. (201 Id) says that this should correspond to 
an SZ signal of T 500 = (51.4 + 4.6) x lO^^arcmin^, much higher 
than our actual tSZ measurement, a conclusion that is indepen¬ 
dent of the mass bias. 

The cluster halo gas scaling relations appear therefore to be 
violated in the QSO environment. This could be due to feedback 
from the QSO SMBH. It is likely that X-ray emission from the 
QSO itself also affects our X-ray measurements. Moreover, we 
are extrapolating the relation to redshifts much higher than those 
used to establish its validity. The redshift evolution of these rela¬ 
tions is in fact not well known, and we have adopted self-similar 
evolution as a reasonable assumption given the lack of empirical 
constraints. Measurements of these relations on non-QSO sys¬ 
tems at the same redshifts would provide a valuable comparison 
sample to judge if the QSO’s activity is in fact responsible for 
violating the scaling relations. 

Einally, we note that interpretation of the measured tSZ sig¬ 
nal is also clouded by possible contribution from gas in struc¬ 
tures correlated with the QSO halos and falling with the rela¬ 
tively large effective beam of Planck, as suggested by Cen & 
Safarzadeh (2015). We leave careful modeling of this effect to 
future work. 

7. Conclusion 

Planck"s coverage of the SDSS survey area gives us the opportu¬ 
nity for in-depth study of quasar environments through submil¬ 
limeter observations of optically selected QSOs. This spectral 
band provides valuable information on star formation and en¬ 
ergy production by the QSO central engines. Eiltering Planck"^ 
seven spectral bands from 70-857 GHz with multi-component 
MMEs, we have measured dust, synchrotron and tSZ signals 
from the large BOSS QSO sample spanning redshifts out to 
0.1 < z < 5. While measurements of individual QSOs are well 
below Planck’s sensitivity, the MMEs enable us to determine the 
population’s mean properties as a function of redshift. 
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Dust emission dominates the average QSO spectrum, and its 
evolution with redshift suggests that star formation in QSO en¬ 
vironments follows the general cosmic star formation rate with 
a peak at z ~ 2. We did not detect a tSZ signal from QSOs at 
z < 2.5, but did hnd a signal at higher redshifts. The signal is 
clearly seen out to z ~ 4, the highest redshift to which the tSZ 
signal has been measured to date. A non-negligible fraction of 
the gas’s thermal energy could be supplied by QSO feedback, 
but interpretation is difficult because of modeling uncertainties 
associated with halo gas scaling relations. Finally, we observe 
synchrotron emission from a radio-loud sub-sample at rest-frame 
frequencies v > 100 GHz and whose luminosity rises monotoni- 
cally with redshift. 

These results are based on a large sample of nearly 300 000 
BOSS spectroscopic QSOs falling outside of the Planck mask. 
Our findings present a fresh and representative view of the aver¬ 
age properties characterizing the QSO population. 
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